Ground-state energy and frustration of the Sherrington-Kirkpatrick model and 

related models 
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Exact ground states are calculated for the Sherrington-Kirkpatrick (SK) spin-glass containing up 
to A'^ = 90 spins. A ground-state energy per spin eff^ = —0.7637 ± 0.0004 is found from the 
dependence of the misfit parameter, which is a measure of frustration of the system. The results are 
compared with those of two related models, which can be introduced by replacing all interactions of 
the SK model by ferromagnetic or antiferromagnetic ones of the same strengths. A parameter x is 
introduced, which describes the fraction of antiferromagnetic interactions in these types of models. 
From the x dependence of finite models it is concluded, that the the SK model {x = 0.5) assigns a 
transition between a ferromagnetic state {x < 0.5) to a spin-glass state {x > 0.5). 

PACS numbers: 02.10.Ox, 05.50.-|-q, 75.10.Hk, 75.10.Nr, 75.40.Mg, 75.50.Lk 



The Sherrington-Kirkpatrick (SK) model was intro- 
duced to describe the spin-glass phase, which is relevant 
to understand the physics of glassy matter. Its impor- 
tance arises from the possibility to get analytical results 
(see e.g. 0), but it is also an object of mathematical 
interest . Currently, the question is controversially de- 
bated, whether spin glasses are 'still complex' after some 
decades of intensive studies (Q, 0, |^ and references 
therein) . Another problem concerns if the results of the 
SK model are valid also for more realistic short-range spin 
glasses 0. In this context, the question of the finite-size 
scaling is widely discussed. 

The problem of finding the exact ground state of most 
of the Ising spin-glass models belongs to the class of NP- 
hard problems, i.e. no algorithm is known which finds the 
optimum in polynomial time @. Therefore, in the past, 
investigations of the size dependence of the ground states 
often were based on approximations. Very recently, at- 
tention has been paid to the problem of fluctuations of 
spin-glass ground-state energies 0, C3- 

In this letter a numerical investigation of exact ground 
states up to A'^ = 90 spins is presented, which are ob- 
tained with the branch-and-bound algorithm of discrete 
nonlinear optimization see also p^ . 

Starting from the Hamiltonian 



n = 



l<i<j<N 



(1) 



Gij, where the 



the SK model is described using 

independent identically distributed random numbers Gjj 
are chosen from a Gaussian distribution with zero means 
and variance one. The chosen scaling in (1) secures 
that the spin-glass ground-state energy is an extensive 
quantity. With Jij = \Gij\ a fully connected ferro- 
magnetic system with the ground-state energy per spin 
e}f = {N - l)/(27rA^)i/2 obtained. Because e}f scales 



like N^/^, the ferromagnetic model should be considered 
in the following as a limiting case of the SK model irre- 
spective of its physical meaning. 

Another limiting case is the fully connected antiferro- 



magnetic (afm) system: J^j 



\Gij\. A further interpo- 



lation of the described models can be obtained consider- 
ing a distribution of Jij according 

P(Jy) = xS{J,, + |Gy|) + (1 - x)SiJ,j - |Gy|). (2) 

Starting from the SK model {x — ^) with (2) the be- 
havior of the system can be investigated, when randomly 
selected ferromagnetic interactions are replaced by anti- 
ferromagnetic ones of the same strengths or vice versa. 

The main feature of the ground state is frustration. It 
is shown in that for Ising spin glasses a 

misfit parameter can be introduced as a measure for 
frustration: 



A^o = ^ (1 - eo/elf) 



(3) 
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It takes into account that starting from an 'ideal' un- 
frustrated system with the ground-state energy Eq''' each 
unsatisfied bond leads to an increase of the ground-state 
energy double its strength 17] . Consequently, e}f in (3) 
is identical to the value for the ferromagnet case {x = 0) 
for all models (2). 

The misfit parameter (3) describes the fraction of each 
bond of the system, which is on average not satisfied in 
the ground state. For the antiferromagnetic triangular 
lattice, for example, the value is hq = i, because one of 
three bonds of equal strength cannot be satisfied. The 
same parameter was used by Stein et al. p3 to charac- 
terize the ground-state energy of ±J models in depen- 
dence on the concentration of antiferromagnetic bonds. 
In a similar way, a related parameter 'ground-state en- 
ergy per bond' is used by Vogel et al. dOj. 

Obviously, the maximum of /lo is ^ for highly frus- 
trated systems (e.g. for high-dimensional hypcrcubic and 
fee fully frustrated ±J systems 0, 0|, 0|). Because 
for the SK model cq = 0(1) and e}f = 0{N^/'^) for 
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FIG. 1: Mean misfit parameter (/Iq as a fmiction of N~^'^ 
for the SK model (x) and tlie fully connected antiferromag- 
netic system (+); error bars are smaller than symbol sizes. 
The lines represent the fits on the basis of formula (3), see 
text. The data points for N = 99 (120,000 samples) and 199 
(64,000 samples) (*) are estimated from Palassini's results 
obtained by a hybrid genetic algorithm jlC|. 



A'' ^ 00, it belongs also to the class of systems with 
maximum occurring frustration (/xq = 5)- (Erroneously, 
the correct N dependence of e^^ was not taken into ac- 
count in [31 , so that the discussion about the misfit of 
the SK model is wrong in that paper.) 

Exact ground-state energy and misfit parameter were 
determined by using the branch-and-bound algorithm for 
the SK model and the fully connected afm model, the 
sizes studied are iV = 19 to 90 and 56, respectively. The 
numbers of samples range from 210,000 {N = 19) to 34 
{N = 90, SK) and to 99 {N = 56, afm). The mean misfit 

parameters (/ip^'') are plotted in Fig. 1 as a function 
of and fitted by (3). For this purpose e}f in (3) 

is expanded in ascending powers of N^^^"^. Assuming 
that the N dependence of the mean ground-state energy 
per spin follows {ef^) = eg° 6A^"" with w = |, the 
fit resuhs in eg° = -0.7637 ± 0.0004, which is close to 
the analytical resuh (e^-^-^ = -0.76321 22]). A fit with 
w = 0.671 shifts even closer to e§^^ . This tendency 
was also found in 10]. A determination of the finite- 
size scaling of the standard deviation a = (((ep^'')^) — 

(gq^'')^)^^^ is restricted by the small number of samples 
for larger N. It results in p ~ 0.72 assuming a oc N^^, 
which is close to Palassini's result (p ~ |) is 
pointed out in that the true values of p is probably 
slightly larger than |. 

Much less is known about the related afm systems, 
which are even higher frustrated already for smaller N. 
Preliminary results for (gq^') vs. N seem to exclude 



FIG. 2: Misfit parameter fio as function of the concen- 
tration X of antiferromagnetic bonds for some systems with 
iV = 49(-^-),56(x),81(★),90(□), 102(A). The iV ^ 00 behav- 
ior is suggested by the straight line (Equ. (4)). 



the validity of the same uj value as it is proposed for 
the SK model. Instead, u — 2 is chosen leading to 
6o°afm = —0.474. The estimation of the N dependence 
of the standard deviation results in p^f,,, = 0.74 ± 0.03. 

The misfit parameter depending on the fraction of 
antiferromagnetic bonds x of the model (2) is shown 
in Fig. 2 for some finite systems. For x <C 1 the 
ground-state energy increases linearly with x according 
to eo(a;) = e'^Q{l — 2x), i.e. Ho{x) = x, as long as the 
ground state remains ferromagnetically ordered. Ob- 
viously, this behavior persists also with increasing x, 
whereas fiQ is stabilized for x > 0.5. (A relatively flat 
minimum for the finite systems can be understood keep- 
ing in mind that the basic element of the system is the 
triangle configuration, which has the tendency to reduce 
the frustration, when it is built of two antiferromagnetic 
and one ferromagnetic bonds. This situation arises rela- 
tively often at x = |.) Otherwise, it can be recognized 
from the results presented in Fig. 1 that lJ,o{x > ^) = ^ 
in the thermodynamic limit N ^ 00. So the numerical 
results for system sizes up to = 102 suggest 

/ ^ f ^ r < 0.5 

M^)-[q^ for. (4) 

and imply that starting from the SK model {x — im- 
mediately a ferromagnetic ground state appears, when 
antiferromagnetic interactions partly are replaced by fer- 
romagnetic ones. 

In summary, exact numerical data for the ground-state 
energy and its finite-size scaling are presented, which are 
in agreement with recent results of Palassini 10] obtained 
using hybrid genetic algorithm. They also provide the 
T — > limit of the energy Cquench of inherent structures^ 
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i.e. the mean energy of the minima accessible from equi- 
librium configurations |2^, Moreover, exact results 
are important to control the efficiency of approximated 
algorithms, which can be applied for larger systems f25|. 

The SK model belongs to the class of maximally frus- 
trated systems, which have a misfit parameter /io = 5- It 
is embedded as a special case into related models, which 
are introduced by varying the signs of the interactions 
Jij. The ground state of the SK model switches over 
from a spin-glass state to a ferromagnetic one, when the 



fraction of the ferromagnetic interactions is larger than 
the fraction of antiferromagnetic ones. The latter conjec- 
ture may also enlighten the controversy on the complex- 
ity of the SK model or its vanishing and is challenging 
for further investigations. 
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assistance of Ch. Roeper in numerical calculations is 
gratefully acknowledged. 
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